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Diffusion Length of Thermal Neutrons in Water 

James A. Dejuren and Hyman Rosenwasser 

The thermal-neutron diffusion length in a hydrogenous medium can be simply deter- 
mined by measuring the variation of the thermal-neutron density at large radii from the 
center of a photoneutron source. The necessary requirement is that the relaxation length 
of the source in a medium be less than the diffusion length, so that diffusion dominates at 
large radii. Using a Ra-Be photoneutron source, the diffusion length of thermal neutrons 
was determined by indium-foil measurements at large radii. A value of 2.763 ± 0.015 cm 
was obtained. This result was confirmed in a second experiment with an Sb 124 -Be source. 



1. Introduction 

In connection with the National Bureau of Stand- 
ards program of absolute calibration of a standard- 
neutron source, the diffusion length of thermal 
neutrons in water was measured. A wide spread of 
values has been published, and the present experi- 
ment is compared with some of the earlier results. 

For a high- energy neutron source the thermal- 
neutron density in water at large radii from the source 
decreases in accordance with the relaxation length 
defined below. If a source having only low-energy 
neutrons is employed so that the rate of neutron 
thermalization per unit volume falls off much more 
rapidly with increasing radius than the thermal- 
neutron density from a hypothetical thermal-point 
source, the thermal-neutron density at large radii 
will vary as a function of the diffusion length and 
radius only (as in the case of the thermal source). 
Photoneutron sources have sufficiently low maximum 
neutron energies for this purpose, and two sources 
having different neutron-energy maxima were selected 
for measurements. The general diffusion equation is 



dn/dt=q-a+DV 2 n, 



(1) 



where n is the thermal-neutron density, q is the rate 
of thermal-neutron production per unit volume 
(thermal-neutron source-strength density), and a is 
the absorption rate per unit volume. 

a=nv/\ a ; D=\ tT v/3, 

where X a and \ tT are the thermal -neutron absorption 
and transport mean free paths, respectively, for the 
average Maxwellian velocity, v, and D is the diffusion 
constant. The steady-state equation is 

V 2 n-n/L 2 +q/D=0, with Z 2 =\ a A, r /3, (2) 

where L is the diffusion length. This equation has 
been amply treated by Fermi and others in standard 
texts. 

With water as the medium and at large distances 
from the source, one can assume with sufficient 
accuracy that 

q=Ke- r/b /r 2 , (3) 

where if is a constant, r is the distance from the 
source center, and b is the relaxation length. 



This relation is very nearly obeyed even for a 
source with a broad spectral energy distribution such 
as is generated by Ra=Bc(a,n). At large r, b is 
usually evaluated by measuring resonance activities 
at energies near thermal with Cd-covered foils. 
Then at large r 
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Since the distribution has spherical symmetry, let 
n=u/r. Then 



V 2 n= lfd 2 u\ 
r\dr 2 ) 



and 



d 2 u u K 
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The solution to the equation satisfying the bound- 
ary condition u-^0 as r-> oo is (derivation in appen- 
dix): 



KL -l 

u=nr= 1 rjje L 
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C- —dx + e^\ —dx \, (5) 
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where ai = l/b — l/L and a 2 — 1/6 + 1 /£. The expo- 

/"» CO 

nential integrals (e~ x /x)dx=EI( — ar) may be 

J ar 

evaluated from "Tables of sine, cosine and expo- 
nential integrals," NBS Mathematical Tables MT5 
and 6 (U. S. Government Printing Office, Washing- 
ton, D. C, 1947). 

As r— > oo t the solution may be expressed as 



Cie 



-r/L 



b 2 L 2 k 



,-r/b. 



(6) 



If i>6, the first term will dominate at sufficiently 
large radii (the necessary requirement for evaluation 
ofl). 

Since the second term acts as a perturbation, 
measurements with sources having different values of 
b serve as a check on the effect of the perturbation at 
the radii where measurements are possible. 

For a Ra-Be(7,n) source the maximum neutron 
energy emitted is only 0.7 Mev (compared to 13 Mev 
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Figure 1. Resonance data for the determination of the relaxa- 
tion length. 

a, Ra-Be (7, ft) source, indium foils; b, Sb 124 -Bc source, gold foils. 

for an a,n source), and the slowing-down length is 
short. Beyond 8 cm from the center of the National 
Bureau of Standards photoneutron standard [1] 1 
the indium resonance activity decreases in accordance 
with eq 3, with a relaxation length 6=2.45 cm (fig. 
1, a). 

2. Experimental Arrangement 

The Bureau photoneutron source was placed in the 
center of a 42-in. diameter stainless-steel tank filled 
with distilled water. Four indium foils, 2.0 by 3.4 
cm 2 in area and 76 mg/cm 2 thick, held in light 
aluminum frames, were suspended by wires from a 
metal framework attached to the top of the tank. 
The foils were centered in the horizontal plane 
passing through the source at positions 90 deg apart. 
A cathetometer mounted horizontally above the tank 
made it possible to position the foils accurately to 
within a fraction of a millimeter at the radii em- 
ployed. The source, a 4-cm-diameter Be sphere 
enclosing a 1-curie RaBr2 capsule, was pin-mounted 
on a vertical stainless-steel rod extending from the 
tank floor. 

Periodic measurements of the water temperature 
gave a mean value of 22° C, with maximum devia- 
tions of 2 deg. 

3. Counting Techniques 

The foils were activated to saturation at radii 14, 
16, 18, 20, and 22 cm from the origin with and without 
0.020-in. Cd covers. All four foils were simultane- 

1 Figures in brackets indicate the literature references at the end of this paper. 



ously counted with "thin" mica end-window Geiger 
counters, the outputs of which were fed to a mixing 
circuit, which permitted the foils to be counted as a 
unit by a single scaler. The sides of the foils facing 
the source also faced the counter windows when 
placed in their respective holders inside the lead 
Geiger-counter shields. Uranium standards were 
used to monitor the Geiger counters, and repetitive 
cycles were made of standard, background, and foils, 
with counting intervals of 5, 10, and 20 minutes, 
respectively, during the decay of the foils. 

4. Results 

The data are given in table 1. A correction of 6 
percent must be made for the absorption of 1.44 ev 
In resonance neutrons by the Cd covers. Errors 
given are determined by the consistency of the data 
as well as the counting statistics. The activities 
given are the counting rates of the foils at saturation 
normalized to the same monitor reading. Figure 2 
is a semilog plot of radius times thermal activity as a 

Table 1. Data from NBS standard Ra-Be photoneutron source 



Position 
(r in cm) 



14.0 
16.0 
18.0 
20.0 
22.0 



Activity 
without Cd 

(Counts 
per minute) 



7,657 ±38 

3, 546 ±18 

1,533 ±9 

655 ±10 

292 ±5 



Resonance 
activity X 
1.06(Counts 
per minute) 



190 ±7 
64 ±15 

a 22 ±6 
*8 ±3 

*3 ±1 



Thermal 
activity 

(A th ) 

(Counts 

per minute) 



7,467 ±39 

3, 482 ±23 

1,511 ±11 

647 ±10 

289 ±5 



104, 540 
55, 710 
27, 200 
12, 940 
6,360 



» Calculated. 
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Figure 2. Indium foils in water. 
Ra-Be (y,n) source. 
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function of r. Beyond 16 cm, rA th decreases 
exponentially. 

^4 r l /i i \ „ T 2 r l 

__ (1:t€)=e ___, 

where e is the experimental error of the ratio 

Using the data at 16 and 22 cm as an example, 
, 5571 , 6 

In w ±€== z 



Table 2. Data from Sb 1?4 -Be photoneutron source 



Z=2.765±0.023cm. 



=0.018 



The weighted least-squares fit of the data between 
16 and 22 cm yield Z=2.763± 0.015 cm. 

5. Sb 124 -Be Experiment 

An antimony pellet that was irradiated in the Oak 
Ridge National Laboratory pile for 2 months was 
substituted for the RaBr 2 capsule in a 4-cm-diameter 
beryllium sphere. The 60-day Sb 124 has three 7-rays 
above the Be photoneutron threshold, the principal 
one having an energy of 1.7 Mev and giving rise to 
24-kev neutrons [2]. As the neutron energies are 
considerably lower than those from the NBS standard 
source, moderation occurs closer to the source, and 
the relaxation length is shorter. A value of 2.1 cm 
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Figure 3. Indium foils in water. 
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Position 
(r in cm) 


Activity 
without Cd 

(Counts 
per minute) 


Resonance 
activity X 
1.06 (Counts 
per minute) 


Thermal 
activity 

(A th ) 

(Counts 

per minute) 


rAth 


14.0 
16.0 
18.0 
20.0 


625.8 ±6 
262.3 ±10 
110. 7 ±6 
50. 6 ±5 


12.4 ±3.0 

a 1. 5 ±0. 5 

a 0. 5 ±0. 2 

». 2 ±0.1 


613.4 ±7 
260. 8 ±10 
110. 2 ±6 
50.4 ±5 


8,588 
4,173 
1,984 
1,008 



a Calculated. 

for b was obtained from the Sb-Be source in water, 
the measurements being made with Cd-covered 
indium and gold foils at radii beyond 8 cm from the 
source (fig. l,b). 

Indium foils that had been intercalibrated in the 
same slow neutron flux were simultaneously exposed 
to saturation at different radii from the source. The 
foils were successively removed and counted with the 
same end-window tube. For radii greater than 12 
cm, rAm decreases exponentially. At 14 cm the 
corrected cadmium reading is 1 percent of the thermal 
reading and decreases to 0.35 percent at 20 cm 
(calculated by using 6=2.1). The least-squares fit 
of the data (table 2) gives a value of L equal to 
2.76 ±0.07 (fig. 3), which is in agreement with the 
Ra-Be result, even though the ratio of indium 
resonance to thermal-neutron counting rates is 
considerably lower for the Sb-Be over the range 
covered. 

6. Discussion and Comparison of Results 
For a purely thermal source located at the origin 
n(r)=Ce- rtL lr 



and 



/» CO 

n (r) r*dr 
=6L 2 . 



x 



n (r) r 2 dr 



In analogy, the slowing-down distance, L s , of 
neutrons from a high-energy source is usually esti- 
mated from the mean-square distance_obtained from 
resonance-neutron measurements. r 2 e ,. = 6Z 2 , and 
the approximate relation 



;=rL. + 6L 2 



(7) 



Sb 124 -Be source. 



is used to calculate i 2 . 

Equation (7) treats the resultant thermal-neutron 
distribution as though the high-energy neutrons 
diff used from the source with a diffusion length equal 
to -yfLf+I/. The shape of the radial thermal and 
resonance-neutron density distributions are actually 
Gaussian near__the origin. As a result, the thermal 
and resonance r 2 values give somewhat overestimated 
diffusion lengths, and the value of L obtained by 
this method should represent an upper bound to the 
true figure. 
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Table 3 lists some of the early values obtained with 
RaBe(a,^) sources, using eq 7. 



Table 3. 



Detector 



Rh[3] 

In [4] 

Dy+In [5] 
In [6]..,.._ 
In [7] 



cm 2 

327 

337 

330 

371 

363 



cm 2 

277 

272 

278 

288 

313 



cm 
2.85 
3.3 
2.94 
3.6 
2.88 



cm 
9.4 
9.4 
10.1 
9.4 



ia 



cm 
6.8 
6.7 
6.8 
6.9 
7.2 



Since eq 7 involves a small difference between two 
numbers, a 1 percent error in either r 2 term causes 
about a 3-percent error in L. 

Berthelot, Cohen, and Reel [8] have determined L 
with a Ra-Be(a,?i) source in a water tank by a 
different technique. A thin spherical shell, which 
could be filled alternately with H 2 (or H 2 0+H 2 S0 4 ) 
and a CdS0 4 solution, surrounded the source. The 
difference in the neutron densities outside the sphere 
for the two cases arises from the diffusion of thermal 
neutrons present when the sphere is filled with water. 
Measurements with foils gave a best value of 2.77 ± 
0.04 cm for the thermal-diffusion length. 

Larsson [9] recently has made extensive measure- 
ments of the same type at different CdS0 4 concen- 
trations. An attempt was made to correct for the 
perturbation in the slow-neutron density in the 
immediate vicinity of the sphere and for the attenu- 
ation of high-energy neutrons by the cadmium. He 
obtained a diffusion length of 3.2 5 ±0.1 3 cm. How- 
ever, a different interpretation of the data would give 
a lower value. 

Perhaps the simplest method for measuring L is 
afforded by a nuclear reactor. If a water tank is 
placed on top of a thermal column, the neutron den- 
sity, a few transport mean free paths from the bound- 
ary should decrease as e~ z,L , providing that the di- 
ameter and height of the tank are much greater than 
the diffusion length. Early reactor experiments 
gave a value of about 2.8 cm. Sisk [10] has recently 
obtained a value of 2.67±0.02 cm at 24° C, using the 
Oak Ridge pile and different boundary conditions. 

Probably most of the variation in the values of L 
obtained by the different techniques is due to the 
failure to meet experimentally the rigorous boundary 
conditions required. 

In the present experiment an upper limit to L is 
obtained unless the bracketed term in eq 5 can be de- 
termined. The exponential integrals have been 
evaluated, using Z=2.76, and G was roughly deter- 
mined for the radium source by calibrating the foils 
in a standard flux. A perturbation of less than 5 
percent ((7>1.3) in L should be caused by neglecting 
the exponential integrals. For the Sb-Be source, 
the exponential integrals are much smaller, and neg- 
lecting them should affect the calculation of L less. 
The Sb-Be (y,n) source was fairly weak, and a stronger 
source would permit measurements beyond 20 cm, 
where the exponential integrals would have negligible 
effect. 
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8. Appendix 



Solution of diffusion equation: 



d 2 u u__ 
dr*~~L 2 ~ 



K 



Let p— 



dr 



K 



(p-l/L)(p + l/L)u=-f r e-^ 

__ -1 K _ lh 

U (p-HL)(p + l/L)Dr e 

Solving by partial fractions : 

1 Si 



7 



(p-l/L)(p + l/L) p-l/L l p + l/L 
£i = L/2, £ 2 =-L/2 

u=LK/2D [-^yjj Vfr 1 

Lp + l/L p-l/LJ 

u = Cie-r' L + C 2 e +r < L + KL[2D 



e~ T ' h 
r 



[r e -r(i/b-i/D 
e -rjL d r _ e r/L I 



g-r(l/6+l/L) 



dr 



Let ai =l/6-l/L, a 2 =l/b + l/L. 



The indefinite integrals may be expressed in terms of definite 
exponential integrals 

f e -ar /*oo e - ay poo ^ 

dr=— dy=— — dx- 

J r Jr y Jar x 

To satisfy boundary conditions at infinity C 2 —0, and 
setting d=KLC/2D. 

2D L J«ir X Jazr X J 

To evaluate C, K must be determined. Since less than 1 
percent of the neutrons from either source are captured 
above thermal energies 

/» oo n 

Q = 4x q(r)r 2 dr = 4:Tr\ —Ke~ r ^ b dr, 
Jro Jro 
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where Q is the neutron emission rate of the source. 

The expression q=Ke~ rn /r 2 does not hold down to the 
radius of the source, so a larger value of tq than the source 
radius must be used. For the Ra- Be (7,71) source r$ is esti- 
mated to be between 4.0 and 4.5 cm from previous data. 
Using K=Qe r <>lb/4:Trb > an absolute measurement of n would 
permit a rough determination of C as all other quantities 
are approximately known. 

Exponential integrals for L=2.76: 



Source 


«i 


<*2 


Ra-Be 


0. 0458 
.114 


0. 7705 
.8385 


Sb-Be 





Ra-Be(7,n) source 


r 


/• 00 

I e-'/xdx 
J air 


€ 2t/l 1 er*/xdx 


cm 
16 


0. 351 
.297 
.253 
.217 


0.037 
.030 
.024 
.019 


18 


20 


22... 






Sb-Be(7,7i) source 




14 


0. 0868 
.0625 
.0455 
.0339 


0.0161 
.0113 
.0078 
.0054 


16 


18 


20 





Washington, July 1, 1953. 
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